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We report experimental observations of hexaroll chaos in inclined layer convection. Two separate 
populations of defect complexes characterize the defect turbulent state, one the classical penta- 
hepta defect of hexagonal planform patterns, and the other a short-lived complex capable of self- 
annihilation. By measuring the defect statistics we show that short-lived defect complexes give rise 
to linear defect destruction rates. This observation explains previous results in other defect-turbulent 
pattern forming systems. 
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Systems that are driven out of equilibrium often show 
similar patterns although the underlying processes can 
be quite different 0|. One challenge is to find measures 
that can quantitatively assess the similarity of different 
patterns. A convenient starting point for a classification 
are topological defects, as has long been known in the 
field of fingerprint recognition yK Much of the previous 
work H H H II II EM ES El IH 113 is concerned 
with the particle- like appearance of defects and with their 
chaotic spatiotemporal dynamics, i.e., defect turbulence. 

A key feature of defect turbulence is the reaction dy- 
namics between different types of defects. Gil et al. |3| 
showed that the topological property carried by defects is 
conserved, causing them to behave as charged particles 
in reactions. This leads to a generalized Poisson dis- 
tribution for defect populations, and reaction rates that 
obey the mass action law, where defects are created with 
rates independent of N, the number of defects present, 
and are annihilated at a rate proportional to N 2 . This 
result was confirmed 5] in a thermal convection experi- 
ment, where they also accounted for defects entering and 
leaving through the boundaries. Recently, however, two 
results have been brought forward that appear to contra- 
dict the notion that mass action is followed. The first of 
these is a Swift-Hohenberg type model of rotating ther- 
mal convection yj, and the second is a calculation of 
the complex Ginzburg-Landau equation with noise added 
fl4| . Both produce defect statistics where defect annihi- 
lation rate has a part proportional to the number of de- 
fects present in addition to the expected N 2 dependence. 
To date an explanation of this phenomenon is lacking. 



Here we report experimental results from non- 
Boussinesq inclined layer convection (ILC) where a in- 
clined layer is heated at one side and and cooled at the 
other. It is well known [lfij that a shear flow is set up 
with a flow that is up at the warm side and down at the 
cold. For non-Boussinesq conditions in a horizontal layer 
a hexagonal planform would be observed Q. However, 
here due to the shear flow one roll orientation is preferred 
leading to a pattern that is similar to the hexarolls de- 
scribed for cylindrical rotating convection 0,^3- ^ n our 
case, we observe a defect turbulent state of spatiotempo- 
ral chaos (STC) (Fig. ^) that we term Hexaroll Chaos. 
We show that hexaroll chaos has two separate popula- 
tions of defect complexes, one the classical penta-hepta 
defect of hexagonal planform patterns, and the other a 
short-lived complex capable of self annihilation. By mea- 
suring the defect statistics we show that short-lived defect 
complexes give rise to linear defect destruction rates ex- 
plaining the earlier observations in other pattern forming 
systems. Hexaroll chaos can be described by amplitude 
equations of the form 



d t A x = g x A x - rnA* 2 A* 3 - AiJ^axjlAjl 2 , 
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dt A 3 = a 3 A 3 - 773^1^2 ~ ^3 J2 a v\ A i 
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Here the convection amplitude is a linear superposition 
of the three modes: V{x) = 2j=i ^-i exp[iqjx] + c.c, 
where the vectors q define a hexagonal lattice, q! = 
(^q,-q/2),q 2 = (0,g),q 3 = (-^q,-q/2), and the 
vector x represents position. The convection cell is tilted 
upward along the x = (1,0) direction, breaking the sym- 
metry of the Ai mode, but the other two remain symmet- 
rical, so the following coefficient relations can be inferred: 

Vl = Vh a l2 = a 32; a H = a 33' a 13 = a 31 ( 2 ) 

where a star denotes complex conjugate. The linear 
growth rates are given by: 

<7i = e - Ae + iuj, a 2 = e,a 3 = a{ (3) 

A rj-i A rri 

where e = — — - is the reduced control parameter, AT 
is the temperature difference driving convection, and AT C 
is the value of AT at onset of convection. u> > repre- 
sents the downhill drift of modes one and three. The 
implication of these equations is that convection rolls of 
mode two appear at e = and the other two modes ap- 
pear at e = Ae. We observe Ae = 1.0 x 10 -3 ± 10~ 4 
which agrees well with the value of 9.7 x 10~ 4 obtained 
from numerical solution of the Navier-Stokes equations. 
Dislocations primarily arise in modes one and three in re- 
gions of low amplitude. These can be paired into defects 
of two distinct types. One type is the classical penta- 
hepta defect (PHD), where a dislocation of one topolog- 
ical sign in one mode is bound to one of opposite sign 
in the other mode. Additionally, we observe a new type 
of defect complex where short-lived pairs of dislocations 
of opposite sign are found in the same mode. We term 
these dislocation pairs same-mode complexes (SMCs). 

We use a gas convection apparatus [13 with a work- 
ing fluid of C0 2 at (4.240 x 10 6 ± 10 3 ) Pa regulated to 
±300 Pa and the mean temperature of (31 ±0.05)°C reg- 
ulated to ±5 x 10~ 4o C. We tilt the plane of the fluid 
layer to 9 = 5° from the horizontal. The fluid layer 
is d = 308 ± 2/im high, and the vertical diffusion time 
t v = 0.56±0.01 s. Prandtl number is a = k/v rs 1, where 
k is thermal diffusivity and v is kinematic viscosity. We 
find onset to rolls at AT C = 20.68 ± 0.01°C and onset to 
hexarolls drifting downward with u) — 0.8t v s/rad above 
that at Ae. We choose a region of interest 142<i x 95d, 
with the long edge parallel to the direction of tilt, lim- 
iting spatial inhomogeneity of e to < 10% during the 
experiment. The large value of AT C makes convection 
strongly non-Boussinesq, with non-Boussinesq parame- 
ters (see, e.g., Ref. GjJ of 70 = 0.1653, 71 = -0.2885, 
72 = 0.2053, 73 = 0.0085, and 74 = -0.2491, where 
70 and 71 are the first- and second-order coefficients for 
the dependence of density on temperature, and 72-74 
are the first-order coefficients for the temperature depen- 
dence of viscosity, heat conductivity, and specific heat, 
respectively. 
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FIG. 1: (a) A sample image from e = 0.054. Uphill direction 
is to right, (b) One oblique mode only, produced by applying 
an angular Fourier filter to (a). Contrast-enhanced for clarity. 

We collected data by quasistatically increasing e in 
steps of 10~ 3 to e = 0.054, where we acquired 10 6 
shadowgraph images of the region of interest at 5 
frames /second, for an observation time of 357100t„. (See 
Ref.GJ for details of shadowgraphic technique.) An ex- 
ample frame appears in Fig. ^ Note the variation in 
oblique roll amplitude that gives the appearance of mixed 
rolls and hexagons, and that oblique- mode dislocations 
are restricted to those areas of low amplitude. 

We processed all images by first separating them into 
three individual modes by applying an angular Fourier fil- 
ter. We identified dislocations in each mode using previ- 
ously described techniques || . Dislocations in mode two 
(with stripes parallel to tilt) are rare, and were recorded 
but not processed further. After all dislocations were 
identified for all frames, we assigned each to either a PHD 
or to an SMC. We did this by iterating through all dis- 
locations, and searching within a threshold radius for an 
opposite charge dislocation in either the same mode or 
the other mode. If such a match was found, both disloca- 
tions were marked as being bound, and the search moved 
on to the unbound dislocations. After all dislocations 
were visited, we increased the threshold and repeated the 



3 



0.2 
0.15 



S 0.1 

o 



0.05 




~i 1 1 r 



n 1 r 

PHD 1= 



smc I-.".---.-.; 



n-,-» | 



01 2345678 
Distance / d 

FIG. 2: Histogram of separation distance of dislocations in 
defect complexes 



search over all unbound dislocations. This process was 
continued until a maximum threshold of ^ 12d is reached 
(Fig. |3J) • Typically, ~ 90% of dislocations are accounted 
for in bound pairs, leaving the other 10% unbound at 
the edges. PHDs are roughly four times more prevalent 
than SMCs in all 10 6 images. The next step of processing 
was to track the center points of PHDs and SMCs by the 
dislocation-tracking method of Huepe et al. .12]. This 
allowed us to tally defect creation and destruction rates 
and rates of entering and leaving the frame (Figs.|31&EJ). 
We calculated the rate for each N by dividing the num- 
ber of events by the total number of frames containing N 
defects, where N is the geometric mean of the two types 
of either PHDs or SMCs. We split the data set into 100 
equal parts to calculate standard deviation of rates. 

Let us consider the types of reactions that SMCs and 
PHDs can undergo. This will illuminate later descrip- 
tions of defect creation and destruction rates. The clas- 
sical case of a symmetrical hexagonal pattern (see, e.g., 
Ref.[H) has six types of PHDs: 1+2", l"2+,l+3-, 1"3+, 
2 + 3~ and 2~3+ where, e.g., 1 + is a positive dislocation 
in mode 1. These follow reactions of the forms 



1+2" 
1+2- 



l - 3 + <-► 
1^2+ <-> 



2~3 H 
0, 



(4a) 
(4b) 



which yield nine possible reactions once permutations are 
considered. Note that all destruction reactions (those 
which proceed to the right) involve two reactants, so from 
mass action we expect a rate that depends on the number 
of PHDs squared. Creation reactions follow a linear rate 
law, with a nonzero intercept due to l|4b|l . Now, if we 
consider the asymmetric case of hexaroll chaos, we must 
disallow dislocations in mode two. There are now two 
types of PHDs possible: l + 3 _ and 1 _ 3+. In addition, 
there are two types of SMCs possible: l+l - and 3+3 _ . 
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FIG. 3: PHD Creation, destruction, leaving and entering 
rates, (C,D,L,E, respectively) along with fits. Number (N) 
is the geometric mean of the populations of the two types of 
PHDs, which are approximately equal. Rate means number of 
events for each N, divided by the number of frames containing 
exactly N. Error bars represent standard deviation, calculated 
by subdividing the data set, and are used for weighting fits. 
Please note that the mean rate is overestimated when the 
standard deviation is on the order of the mean because all 
values are positive. 



The reactions are: 



i+r 



3+3" <-> 1+3" 
1+1" <-> 0, 



1"3 H 



3+3- 



0. 



(5a) 
(5b) 
(5c) 



Where (|5a|l is the realignment of two SMCs into two 
PHDs, and i|5b() and (|5c|l are the self-annihilation of 
SMCs. Here we have two sets of rates, those for PHDs 
and those for SMCs. PHDs are created by l|5a[) proceed- 
ing to the right. The rate for this reaction is indepen- 
dent of the number of PHDs because there are no PHDs 
as reactants, and populations of SMCs and PHDs are 
uncorrelated. PHDs are destroyed by the same reaction 
moving in the opposite direction, with a rate proportional 
to the number of PHDs squared. SMCs are created by all 
three reactions going to the left, which gives a constant 
rate law. Finally, SMCs are destroyed by these reactions 
going to the right, giving a linear rate from (|5b|) and 
(|5cf) with a squared term from (|5a|l . In practice, we find 
that this squared term is negligible, so it is set to zero 
when constructing fits. Otherwise, these rate relations 
and those for entering and leaving given by Daniels and 
Bodenschatz || fit the data in Figs. El & El with good 
results. 

Using detailed balance, we recursively define a modi- 
fied Poisson distribution of the probability of finding TV 
defects in the frame at a given time in terms of creation, 
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FIG. 4: SMC creation, destruction, leaving and entering rates 
(C,D,L,E, respectively). Number (N) is the geometric mean 
of the populations of the two types of SMCs 
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FIG. 5: Probability distribution functions (PDFs) for average 
numbers of each type of SMCs and PHDs. Boxes are mea- 
sured, and curves are calculated from fits in Figs. [3]&2]using 
detailed balance. 



destruction, leaving, and entering rates: 

C(j - l)E(j - 1) 



N 



v(n) = p(o) n 



mm 



(6) 



where 'P(O) is supplied as a normalization condition. 
Note that for the case of SMCs where C and E are con- 
stants, and D and L are linear functions of N, © reduces 
to the Poission distribution. For the PHD case, © re- 
duces to the squared Poisson distribution of Gil et al. 
jj| with the finite-size extension of Daniels and Boden- 
schatz [i| (Fig. 0. We conclude that defect statistics 
of hexaroll chaos may be easily understood within the 



existing framework of reactions and rate laws if we intro- 
duce a new type of short-lived defect complexes (SMCs) 
that usually self-annihilate, but that may stabilize by 
binding into long-lived PHDs. These complexes appear 
to arise because of broken symmetry which changes the 
allowed ways of creating and destroying PHDs. We sug- 
gest that this may also be the case in other systems, such 
as the rotating hexagons of Young and Riecke @ where 
their ansatz of mixed linear and quadratic rate laws may 
be separable into two sets of laws, one for SMCs and 
one for PHDs. Similarly, Wang's STC state with noise 
added 0] may contain a population of short-lived, noise- 
induced defects that primarily annihilate with the defect 
they were created with, thus following a linear rate law. 
More detailed study of these systems is called for to de- 
termine if there is a generality that connects symmetry 
breaking with altered rate laws. 
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